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Abstract 

New Broer-Kaup type systems of hydrodynamic equations are derived 
from the derivative reaction-diffusion systems arising in SL(2,R) Kaup- 
Newell hierarchy, represented in the non-Madelung hydrodynamic form. 
A relation with the problem of chiral solitons in quantum potential as a di- 
mensional reduction of 2-1-1 dimensional Chern-Simons theory for anyons 
is shown. By the Hirota bilinear method, soliton solutions are constructed 
and the resonant character of soliton interaction is found. 



1 Introduction 

Recently, a modification of the nonlinear Schrodinger (NLS) equation by a quan- 
tum potential has been studied in several problems arising in low dimensional 
gravity, [T], plasma physics, [5], the capillary wave, and information theory, [5]. 
Subsequently, the influence of this potential on anyons in 2+1 dimensions has 
been studied [J, and the Abelian Chern-Simons gauge field interacting with NLS 
has been represented as a planar Madelung fluid 6 , where the Chern-Simons 
Gauss law has the simple physical meaning of creation of the local vorticity for 
the flow. For the static flow when the velocity of the center-of-mass motion is 
equal to the quantum velocity, the fiuid admits an N- vortex solution. It turns 
out that in this theory the Chern-Simons coupling constant and the quantum 
potential strength are quantized. 

Reduction of problem to 1-1-1 dimensions leads to JNLS and some versions 
of Derivative NLS in quantum potential. Hence the chiral solitons appear as 
solutions of the derivative NLS with quantum potential ^ . The last one 
by the Madelung transform is represented as the derivative Reaction-Diffusion 
(DRD) system, arising in SL(2,R) Kaup-Newell hierarchy, and giving rise to the 
resonant soliton phenomena [J.. 



In the present paper by using new, the non-Madelung representation, we 
formulate the problem in terms of novel hydrodynamic systems of the Broer- 
Kaup type. Then by Hirota's bilinear method we construct chiral solitons for 
the system and show the resonance character of their interaction. 



2 Dimensional reduction of Chern-Simons the- 
ory 

We consider the Chern-Simons gauged Nonlinear Schrodinger model (the Jackiw- 
Pi model) with nonlinear quantum potential term of strength s 0]: 



L = ^e^^^A^a.^A + ^(^i^o^ - ^^0^) - + sVl^lVlV'l + ^(V^^), (1) 

where Z?^ = + ieA^, [ji = 0, 1, 2). Classical equations of motion are 

zZ)o^ + DV + 1^V = s^V, (2) 

m 

diA2 - d2Ai = -^b^P, (3) 

doAj - djAo = ~-ie,k{i>Dk^ - ^PDk^), (j, fc = 1, 2). (4) 

We consider dimensional reduction of this model when all field are indepen- 
dent of X2 space variable, so that 82 — 0- Then in terms of Aq = Ao + eA"^, and 
B = A2, we obtain 

i{do + ieAo)ij + (di + ieAifi^ + V'^ = 5^1^, (5) 

m 

9o^i - ai^o - (6) 



d^B = (7) 

d^B = i- + ieAi)V' - V'lSi - ieA{)^ (8) 

Here and below we skip the tilde sign for At^. The last two equations are 
compatible due to ([5]) and the corresponding continuity equation 

= id^ [^[di + ieAx)^ - il){di - ieAi){p\ (9) 

implies compatibility of equations (O,®. Integrating these equations we 
find B in terms of density Tpip 

B = - [ ipipdx'. (10) 
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From another side, flatness of connection ([SJ implies = ^qi/), Ai — dicj), 
and these potentials can be removed by the gauge transformation, tp = e~'"^\l/. 
As a result we obtain the Schrodinger equation with self-interacting nonlinear 
potential Vlipip), and the quantum potential 

i9o* + 5?* + l^'^' = s?^* (11) 
1^1 ^ ^ 

2.1 Madelung representation and RNLS 

If we substitute the Madelung Ansatz — yfpe~^^ to wave function in (fTTj) 
then we get the coupled system 

ao5-(9i5)2 + r(p) + (i-s)^ = o (12) 

9op-5i(2p9i^) = (13) 
For velocity field v = —2diS it implies the hydrodynamical system 

dov + vdiv^2di(v'{p) + {l-s)^^^ (14) 

dop + di{pv)=0 (15) 

which is the Madelung fiuid representation of (ITT|) . 

First we consider the under-critical case, when the strength of the quantum 
potential s < 1. Then introducing rescaled the time and the phase t — t\/l — s, 
S = -y^^, for new wave function ^ = y/pe~^^ we get 

V' ~ 

ido^ + df'^i + - — = (16) 

In the over-critical case when s > 1, we can't reduce (|lip to (IT5|) . However, 
if we introduce pair of real functions 

e^+\x,t)^^e^, e(^)(x,i)-Vpe-^ (17) 

instead of one complex wave function, then we get the time reversal pair of 
reaction-diffusion equations 

55e(+) + 52e(+) - ^e(+' = (18) 

- 9ge(-) + d^eS-^ - j^eS-^ = (19) 

where f = ^^^~T, 5*= -^^f^, p^e^+^e^-), V ^ V{e^+'> e^-^). 

If interaction between material particles is the delta function pair form then 
potential V{p) — gp^ /2 and equation ([TT|) becomes 
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ido^ + df^ + g\^\H = s^^* (20) 

We called this equation the Resonant Nonlinear Schrodinger equation (RNLS) . 
It appears in the study of low-dimensional gravity model on a line, the Jackiw- 
Teitelboim model JLj , and in description of cold plasma . For under-critical 
case s < 1 it reduces to the standard NLS equation (|16p , and is integrable model. 
For the over-critical case s > 1 it reduces to the couple of cubic reaction-diffusion 
equations (HH), ((111), which is also integrable as SL(2,R) NLS from the AKNS 
hierarchy. In the last case new resonance phenomena for envelope solitons take 
place [1]. In the next section we will discuss another reduction of Chern-Simons 
theory with dynamical field B, and show that in this case resonant versions of 
the JNLS and DNLS equations appear. 



3 Dynamical BF theory 

To do the gauge field component B in Section 2 to be dynamical, following [11] 
we introduce the corresponding kinetic term so that 



(21) 



L = nB{doAi - diAo) +edoBdiB + ^{ijj{do + ieAo)^P 

-iPido - ieAo){j) - {di - teAi)ii;{di + ieAi)i^ + sdi\i;\di\^P\ + ^(V^V), 
Then equations of motion are 

i{do + ieAn)i: + {di + leA^f^ + V'^: = (^2) 
9oAi - diAo = —dod^B (23) 

K 

d,B = -p (24) 
doB = --J (25) 

K 

where the particle and the momentum density are 

p^{pi}, J^-[{ij{dx+ieAi)i}-%l;{di-ieAi)i}\^j + 2eAip (26) 
i 

and J = —ilipdiip — ipdiip]. Equation (j22p implies the conservation law 

dop + diJ^O (27) 

This conservation law is the compatibility condition for the system p4| . (l25|) . 
and allows us to write 
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dodiB — —[adop — (1 — a)diJ] 



(28) 



where a is an arbitrary real constant. Substituting (1281) to (j23p and combining 
terms under the same derivatives we have 



The system 



do[Ai-—apj-diiAo^^{l~a)j\=0 (29) 
5]) is invariant under the local U{1) gauge transformations 

(30) 



Then solving we have 

26*6 ^ , 296,^ , ^ „ 
Ai = -^ap + di0, ^0 = ^(1 - a) J + daq 

and for the gauge invariant field ^ ~ e"^"^^ it gives 



(31) 



We' 



dm 



(32) 



(33) 



where J = j + ^^ap^, j = - p = Finally we have 

i^t + + [(2a + 1)1*12^-, + (2a - 1)^-2^,] 

+4^a(a - 2)|^|4^ + V'^ - 

where partial differentiation notations are evident. The remaining gauge trans- 
formation for this equation is just the global U{1) transformation: ^I* — )■ e''^^', 
A = const. 



3.1 Reductions of general RDNLS 



The behavior of equation (I33p depends on value of parameter s. If we replace 
^ = e^~^^ then we have couple of equations 



~2 S xC 



29e' 



AaRxc'^ = 0, 



St-Sl + il-.s){Rxx + Rl)+ ^ 
determining the Madelung fluid representation 



a{a - 2)e^''- + y' = 0, 



Pt + {pv + ^^2ap')x = 0, (34) 



Vt + VVx 



(l-s) 



^Pxx 
V~P 



29e' 
-^pv 



29e' 



a{a - 2)p2 + V' 



(35) 
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For s < 1 for redefined variables t^/l~~s = t, S/\/l — s = S, ^ = e , we 



liave 



(2a + f)|*|2^', + (2a-f)'I'2^r, 



(36) 



V 



^' = 



Similar to the Chern-Sinions 2+f dimensional case [4] we have effective 
result of the quantum potential in the rescaling of the statistical parametr 

— > K^^/l — s, but in contrast no quantization of this parameter now appears. 
Transformation between wave functions has nonlinear form 



(37) 



For s > 1 it is impossible to reduce the system to the Schrodinger type 
form. However for redefined parameters t\/s — 1 = t, S/ \/s — 1 = S and two 



real functions = e^+'^, E 



^ we get 



TEf + i?± T 
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{2a + l)E+E~E^ + (2a ~ l)E^^E^ 



4^a(a - 2){E^E-rE^ = 



(38) 



3.2 Gauge transformation 

We notice that in the gauge potential representation ([31]) , the gauge function 
(j) = depends on a: 



Comparison with the case a = 
gives relations 



(39) 



(40) 



(41) 



Compatibility of this system is ensured by the continuity equation (j27p . 
Then corresponding gauge transformed wave functions 'i''-"^ and 



(42) 



are related by 
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Integrating (PT|) and substituting to we have gauge transformation be- 
tween equations ([33| and the same equation with a — 0: 



= exp J%dx'^ (44) 

From this relation we can connect two samples of equation p3p with different 
constants a and /3 

= exp (^-«^(a - (3) J" pdx'^ (45) 
Indeed one can check easily from 



= $(/3)^(/3)^ «'(")(9i + ii/ap)^'("' = *('^)(ai + iz/ap)*('^) (46) 
that = p(^), J(") = J(^), and 

{di + iz.ap)*(") = e"*"("-'^^ Z"' (ai + iiyl3p)¥'^'^ (47) 

(9o + - a)J)¥"'> = e~*"(""'" /"^ (So + i:^(l - (3)J)¥^'> (48) 
where z/ = ^H- 

The gauge transformation P5|) for the Madelung representation implies 

„ s-C^) = ;,(q, - /3) y pdx' + 2?™, = i?^'') (49) 

For s < 1 it gives U{1) gauge transformation for ([55)1 in the form 



^ ^(/^)e-^^("-'^)/^'"'-'e-^7& (50) 

The last multiplier can be absorbed by the global phase transformation on 'J. 

For s > 1 the above U{1) gauge transformation give rise to the local 5*0(1, 1) 
scale transformation (the Weyl transformation) for equation 



4 Integrable DRD Systems 

It was shown above that the one dimensional problem of anyons in quantum 
potential with a specific form of the three-body interaction, can be reduced to 
the general resonant DNLS equation. 



= f pdx' 
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4.1 General Resonant DNLS 

This equation 

+ + iD [{2a + 1)1*1'*,, + (2a - 1)*'*^ 



+4^'("- 2)("- 2)1*1 * = ^W* ^^^^ 



is integrable for any values of parameter a. 



4.2 The Resonant Case 



For special case s > 1, by the Madelung transformation * — e^^*'^ and intro- 
duction of two new real functions = e^~^^ , E~ = e^~^ we get the general 
DRD system 



(2a + l)E+E-E^ + (2a - l)E^^E^ 



^)(« - UE^E-fE^ = 0, 



where d is the statistical parameter. 

This system has particular reductions 
1. DRD-I, (a = 3/2) 



2. DRD-II, (a = 1/2) 



(53) 



-E++E+,-2u{E+E-E+)^ = 0, (54) 
+E^ +E-, + 2u{E+E-E-), = (55) 



-Et + Et,-2vE+E-Et = 0, (56) 
+E- +E-^ + 2vE+E-E- = (57) 

3. DRD-III, (a = -1/2) 

- E+ + E+^ + 2iyE+^E- -2iy^{E+E-)^E+ = 0, (58) 
+E- +E-^-2uE-^E+ -2u^{E+E-fE- = (59) 

4. JRD, (a = 0) 

T i;*^ + - ^ [E^E- - E+E-] E^ 

— —{E+E-fE^ = (60) 
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5 Resonant Hydrodynamic Systems 

To find hydrodynamic form of the above equations we introduce velocity vari- 
ables according to the Cole-Hopf transformation 

v+ = {lnE+%, v- = {lnE-% (61) 

and density 

p = E+E-. (62) 

Then by identity 

Px = pv^ + pv~ , (63) 

we can rewrite the DRD system in a closed form for only one of the couples of 
hydrodynamic variables {p,v~^) or {p,v~). 

5.1 Hydrodynamic Form for DRD-I 

For DRD-I case it gives the new hydrodynamic system 

Pt + Pxx = [2pv+ -3iyp%. (64) 

5.2 Hydrodynamic Form for DRD-II 

For DRD-II case first we get the coupled heat equation with transport 

- E+ + i?+ - 2vpE+ = 0, 
Pt + Pxx = i2p{lnE'^)^ - vp'^)x- (65) 
Then the hydrodynamic form for this system is 

4 = [vt + {v+f -2upv%, 

Pt + Pxx = i^pv"^ - vp\. (66) 

5.3 Hydrodynamic Form for DRD-III 

For DRD-III case it gives the new hydrodynamic system 

4 = [4 + + Mp - - '^^^p\, 

Pt + Pxx = [2pv^ + vp%. (67) 

5.4 Hydrodynamic Form for JRD 

For JRD case it gives the new hydrodynamic system 

vt = K + {v^? - ^i?pv+ - Px) - \^'pX, 

Pt + Pxx = [2pw+]x- (68) 
In all above cases for v~ we have the system with replaced t — )• —t, u ^ —u. 
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5.5 Generic Case 

For the generic case of arbitrary a firstly we have the system 

- E+ + E+, - v[2pEt + {2a - l)p,E+] 
-v^{a~]^){a~^-)p'E+ = 0, 
pt = [2p(ln_B+)a; - px- 2vap'^]x- 
It gives the new hydrodynamic system 

= ["t + {v^f - J^(2pz'+ + {2a - l)px) 
-u'^{a- '^)p\, 
Pt + Pxx = [2/5W+ - 2vap^]x. 



(69) 



(70) 



6 RNLS and Broer-Kaup system 

The RNLS for s > 1 can be transformed to the reaction-diffusion system 

R+ = R+^ + 2vR+R-R+, 



— i?( — Rxx ^" 2vR^ R R 



(71) 
(72) 



By substitution w+ = {\\\E^)x, p = E^E , it can be transformed to the the 
hydrodynamic form as the Broer-Kaup system, [5], |10) . 



= {^t + + ^I'Px, 

Pt + Pxx = {2pv^)x. 

If = {\nE^ )x, p = E^E^ , then we have 

- «r = {vx + {'"^f)x + 2i^px, 

-Pt + Pxx = (2pw~)x- 



(73) 



(74) 



7 Relation with Broer-Kaup System 

Given E^{x, t), E~'{x, t) satisfying general DRD system (|53| . then real functions 



R-" 



R- 



E- + (a - -)vE+E-E^ 



Sa+^)v E+E- 



(75) 



or 
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R+ = 



R- 



-E+ + {a- -)yE+E-E^ 



E+E- 



E-e 



E+E- 



satisfy the reaction-diffusion (RD) system 



R+ = R+^ + 2iyR+R-R+, 
-R7 = R-^,+2vR+R-R- . 



From this fact we can get next result. 

If and pE satisfy (f70| then and pR determined by 



(76) 



(77) 
(78) 



PR 



v% - {a+-)vpE, 
{pe)x - PeVe + (" - h'^PE' 



(79) 
(80) 



is solution of the Broer-Kaup system ([75)). For and pE satisfying the analog 
of system ((70)) . 



PR 



= VE + {a+^)l^PE + [HvE + {(^-\))>^PE] 

= PeVe + (a - ^}'^Pe 



is solution of ((74)) . 

Similar way we can get result. 

If and pe satisfy (1701) then w^jj and p^ determined by 



(81) 
(82) 



PR 



B - + ^^I'PE + [H-^E + ^))vPe]x, 
'PEV^ + {a- hvpE 



(83) 
(84) 



is solution of the Broer-Kaup system ((73| . For and pE satisfying the analog 
of system ((70)) . 



PR 



Vj^ + (a+ -)i^pB, 

-{pe)x + Pe-Ve + (a - \)vp\ 



(85) 
(86) 



is solution of ((71)) . 
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8 Backlund Transformation 

When p = 0, both systems (|70l) and ((73|) reduce to the Burgers equation. Then 
the above Miura type transformations reduce to the auto-Backlund transforma- 
tions 

^'j? = + (Ini'B)x, = + (Inw;;;;):^ (87) 
for the Burgers and anti-Burgers equations correspondingly. 



9 Classical Bousinesque Systems 

If in ([75)1 we change variables 

p+ + 2vp, (88) 

then we get the classical Bousinesque system 

vt = {{v+f+P+).. (89) 
pt = {vl, + 2p+v+),. (90) 

Similar way in (1741) by variable change 

p^ + 2vp (91) 

we get 

-V = {{v-f+p-),, (92) 
-Pt = {v-, + 2p-v-),. (93) 



10 Bilinear Form and Solitons 

By substitution — / to we have bilinear representation 

{TDi + Dl)ig^ ■f^)^0, (94) 

Dlif+-n + lD,i9^-g-) = 0, (95) 

■ f-)+ag+g- ^0, (96) 

where a = ^ ( DRD-II case), or a = — ^ (DRD-III case). We note that only 
in these two cases the Hirota substitution has simple bilinear form. Then for 
solution of the hydrodynamics systems (|66l) and (|67|) we have 



= {lnE+), = p-^, (97) 
p = E+E- = {In ^),. (98) 
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Bilinearization for arbitrary a can be derived by the gauge transformation, so 
that ^ 

E+ = ^ ,E- ^ ^ (99) 

(/+)^+"(/-)^"" (/+)^-"(/-)^+" 

It imphes next substitution for equation (|70p 

v+ = {lnE+), ^ ^ _ (i + «)|^ _ (i _ a)^, (100) 

f+ 

p^E+E- ^ {In -j^)^. (101) 

10.1 One Soliton Solution 

For one sohton solution we have 

g± = e''? , ^1 + e"*" 6"^ , (102) 

where e"*" — TttttTm?' Vi — kfx ± [k^Yt + ri^^'^\ For regularity of this 
solution we choose conditions /cj~ > and k'^ < 0, then — u < k < v, where 
k — kf + k^ , V — k^ — k^, —kxQ = rj^^"^ + rj^ '"^ + 4>fi- Then velocity is 
positive w > 0, so that our dissipaton is chiral. For the density we have soliton 
solution 

p^E+E-= : (103) 

V — fc^ cosh k{x — vt — xq) + V 

where 2xo = x^ + Xq , and for velocity field 

v+ = ^ ^1 ' r , (104) 

1 + e-^iie''! 

the kink solution 

, V k k 
^ - -teLuh-{x - vt - xo). (105) 

10.2 Integrals of Motion 

The particle number, momentum and energy integrals are given respectively 

f°° 1 f+ 

/ pdx^—\n^\^^ (106) 

J — oo ^ / 

r°° 1 

P=- pv+dx = - Hf+r U-oo (107) 

J— oo 
poo 

E =- [p{v+)'^ ~ Pa:V+ - l^p^V+]dx. (108) 



13 



Then substituting for one soliton solution we find 

1/ V — \k\ V 2v 

The mass of soliton M = \k\/{yv) in terms of particle number becomes M = 
i tanh ^^ , and for the momentum and the energy we have non-relativistic free 

particle form P = Mw, E = 

For the process of fusion or fission of two solitons then the next conditions 
should be valid 

Ar = iVi+iV2, p = p^+P2^ E = Ei+E2 (110) 
Using (|109p after some algebraic manipulations we get the resonance condition 

1^1-^21 = |fcl| + |fc2| (111) 

where Va ^ K ^ K K = K + a = 1,2. 

10.3 Two Soliton Solution 

For two soliton solution we have 

c,± = e"? + e"? + a±e''^+"2"+"? + a±e"^ +"1"+"^ , (112) 
2 

/± = 1 + ^ e^£e''-^+''7 + /3±e"i'+')r+')2++')2" , (113) 

where rif = kfx ± {kf ft + 77*, fcj}'" = (fcf + fc™) and 

+ _ , 1 ^^(fcf - ± _ I kfjk^ - kf)"^ 

' - ^2 (fc+-)^(4^)^' ^ ~ 2 (fc+-)2(4^)^' ^'''^ 



(fel~ fe2')^(^l ^2 )^ ^ 
4(^11 ^12 ^^^21 ^22 



kfkf, (115) 



^2(4-)2' 2(fc+-)^' 2(4-)2- ^^^^^ 

By regularity we have k'l < 0, k^' > in the Case 1, and k^ > 0, k~' < in 
the Case 2. Then solving the resonance condition (jllip we find that for every 
solution of this algebraic equation, the coefficient /3 vanishes or becomes infinite. 
In both cases two soliton solution reduces to the one soliton solution. Hence the 
solution describes a collision of two solitons propagating in the same direction 
and at some value of parameters creating the resonance states (see Fig.l and 
Fig.2). 
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Figure 1: 3D plot of typical soliton resonant state with one soliton resonance 

11 Conclusions 

The problem of chiral solitons in quantum potential, as a reduction of 2+1 
dimensional Chern-Simons theory, was formulated in terms of family of inte- 
grable derivative NLS equations by the Madelung fluid representation. By us- 
ing new, non-Madelung fluid representation we constructed integrable family of 
hydrodynamical systems of the Kaup-Broer type. By bilinear method we found 
resonance character of corresponding chiral soliton mutual interaction. 
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